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SUMMARY 


The  relative  amplitude  method  of  determining  earthquake  source 
orientations  is  applied  to  short  period  teleseismic  array  observations  of  the  1976 
Uzbekistan  earthquake  sequence.  Many  of  these  events  are  shown  to  be 
consistent  with  a  double  couple  source,  and  the  method  provides  well  determined 
orientations  for  most  events  with  magnitude  m^  between  4.0  and  5.6. 
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1.  THE  1976  UZBEKISTAN  SEQUENCE 

This  report  concerns  a  sequence  of  earthquakes  near  the  town  of 
Gazli,  in  Uzbekistan,  USSR  (figure  1),  which  continued  for  several  months  from  8 
April  1976.  The  sequence  included  three  shocks  with  body  wave  magnitude  (m. ) 
above  6,  which  are  discussed  by  Pletnev  et  al.  (1),  and  many  large  aftershocks,  at 
least  sixty  of  which  were  recorded  at  teleseismic  distances.  A  list  of  all  those 
earthquakes  in  this  sequence  which  were  reported  by  the  USGS  National 
Earthquake  Information  Service  (NEIS),  is  given  in  table  1.  Each  event  is  given  a 
reference  number  for  the  purpose  of  this  report.  Table  1  also  shows  their 
magnitudes  as  computed  by  NEIS.  For  the  time  period  of  interest,  there  were  six 
additional  earthquakes  reported  in  the  Bulletin  of  the  International  Seismological 
Centre  (ISC)  -  these  are  shown  at  the  end  of  table  1.  Figure  2(a)  shows  the  ISC 
hypocentres  plotted  on  a  map  of  the  area.  All  but  three  of  these  lie  in  an  area  of 
60  km  square,  but  no  spatial  trend  is  apparent.  In  section  9  this  question  will  be 
re-examined  using  epicentre  relocation  with  the  aid  of  focal  depths  determined 
from  pP  in  this  report. 

This  sequence  of  earthquakes  is  of  considerable  geophysical  interest 
in  that  it  occurred  far  from  any  major  tectonic  zone  in  an  area  of  low  seismic 
activity  (2),  with  no  previous  teleseismically  recorded  earthquakes. 

The  sequence  is  also  of  interest  in  the  context  of  earthquake/explo¬ 
sion  discrimination;  the  wide  range  of  magnitudes  offers  an  ideal  opportunity  to 
discover  how  well  focal  mechanisms  can  be  determined  as  a  function  of  m,  for  a 
large  number  of  similar  earthquakes  in  a  structurally  simple  area.  Our  ability  to 
compute  focal  mechanisms  (ie,  "fault  plane  solutions")  from  teleseismic 
observations  has  been  extended  to  smaller  magnitudes  by  the  relative  amplitude 
method,  and  the  success  of  this  method  introduces  the  possibility  that  some 
events  may  be  positively  identified  as  earthquakes  if  observations  are  found  to  be 
compatible  with  a  widely  accepted  earthquake  source  model.  Such  a  positive 
identification  could  become  particularly  valuable  at  smaller  magnitudes  where 
alternative  data  (for  example,  long  period  surface  waves)  may  not  be  available. 

2.  PURPOSE  OF  THIS  REPORT 


In  this  report,  the  relative  amplitude  method  of  Pearce  (3-5)  is 
applied  to  short  period  seismograms  of  the  1976  Uzbekistan  earthquake  sequence 
as  recorded  by  the  four  short  period  medium  aperture  arrays  at  Eskdalemuir, 
Scotland  (EKA),  Gauribidanur,  India  (GBA),  Warramunga,  Australia  (WRA)  and 
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Yellowknife,  Canada  (YKA).  Phased  array  sums  are  used  throughout  in  order  to 
exploit  the  improvement  in  signal  to  noise  ratio  which  array  records  provide;  we 
expect  a  maximum  signal  to  noise  ratio  improvement  of  /n  where  n  is  the 
number  of  seismometers  included  in  the  sum. 

Application  of  the  relative  amplitude  method  to  this  earthquake 
sequence  should  provide  the  answers  to  two  major  questions.  First,  what  is  the 
practical  capability  of  the  method  to  positively  identify  these  events  as  double 
couple  sources  (ie,  earthquakes)  given  no  more  than  four  short  period  (array) 
stations  at  teleseismic  distances?  Second,  how  well  can  we  constrain  their  focal 
mechanisms  using  the  same  stations?  Of  special  importance  in  relation  to  both 
questions  is  the  magnitude  threshold  below  which  the  capability  is  lost,  or  is  no 
longer  reliable. 

3.  SOME  COMMENTS  ON  THE  RELATIVE  AMPLITUDE  METHOD 


The  relative  amplitude  method  was  introduced  by,  and  is  fully 
described  by,  Pearce  (3, it).  Extensions  of  the  method  and  additional  discussion  is 
contained  in  reference  (5). 

Since  for  shallow  earthquakes  direct  P  and  the  surface  reflections  pP 
and  sP  are  not  well  separated  on  long  period  seismograms,  we  use  short  period 
records,  which  in  any  case  give  us  a  lower  detection  threshold.  It  is  important  to 
note  that  this  relative  amplitude  method  (in  contrast  to  methods  which 
determine  focal  mechanisms  from  P  first  motions)  is  more  suited  to  small  than  to 
large  earthquakes,  since  small  earthquakes  are  more  likely  to  have  direct  and 
surface  reflected  phases  which  are  observed  as  discrete  pulses  at  short  period. 
This  we  expect  because  small  earthquakes  tend  to  have  a  duration  of  source 
radiation  which  is  within  the  short  period  bandwidth  (ie,  less  than  1  s). 

As  used  in  this  report,  the  relative  amplitude  method  assumes  that 
the  source  emits  the  P  and  S  wave  radiations  of  a  double  couple,  and  identifies 
those  ranges,  if  any,  of  double  couple  source  orientations  which  are  compatible 
with  a  series  of  relative  amplitude  observations  made  at  different  stations,  as 
explained  in  references  (3)  and  (4).  Hence,  for  an  ideal  double  couple  source  with 
correct  description  of  ray  takeoff  angles,  source  layer  and  surface  layer 
velocities,  and  with  the  correct  identification  and  measurement  of  each  phase, 
we  expect  the  method  to  yield  at  least  one  acceptable  range  of  source 
orientations  and  we  expect  these  acceptable  range(s)  to  become  smaller  (but 
never  to  disappear)  as  more  observations  are  incorporated  into  the  computation. 
Conversely,  for  a  non-double  couple  source,  again  with  correct  description  of 
structural  and  ray  information,  we  expect  the  acceptable  regions  to  get  smaller 
and  to  disappear  as  the  number  of  observations  is  increased.  Thus,  we  see  that 
the  number  of  acceptable  orientations  will  always  decrease  as  the  number  of 
observations  is  increased,  but  this  tendency  may  mean  either  that  the  source  is  a 
double  couple  which  is  becoming  more  closely  constrained,  or  that  it  is  not  a 
double  couple  and  that  the  inclusion  of  additional  observations  will  eventually 
eliminate  all  compatible  orientations. 

It  may  therefore  appear  that  the  potential  of  the  method  as  a 
discriminant  is  very  limited.  However,  experience  has  shown  that,  if  there  is  any 
error  in  the  representation  of  the  seismograms,  or  if  the  source  radiation 
deviates  from  that  of  a  double  couple,  then  all  orientations  are  eliminated  after 
only  a  small  number  of  observations  is  included  -  say,  two  or  three  seismograms 


of  average  quality  -  whereas  compatibility  persists  to  a  far  greater  degree  when 
the  representation  and  the  source  model  are  basically  correct. 

There  are  two  empirical  results  from  previous  work  which  are 
particularly  relevant.  First,  a  small  deviation  from  ideal  source  radiation  causes 
complete  elimination  of  all  orientations  after  only  a  small  number  of  observa¬ 
tions  is  included.  Second,  it  has  been  found  that,  for  an  ideal  double  couple 
source,  the  region  of  acceptable  orientations  typically  reduces  to  a  small  but 
finite  range  after  only  several  observations  have  been  included  (3,5).  In  cases 
where  all  orientations  are  eliminated,  we  require  to  ascertain  whether  this  has 
arisen  through  an  error  in  the  specification  of  structural  or  seismogram 
parameters,  or  whether  the  source  genuinely  does  not  radiate  as  a  double  couple. 
As  explained  in  reference  (4),  the  relative  amplitude  method  has  provisions  to 
investigate  such  sources  of  anomaly. 

4.  OBSERVED  SEISMOGRAMS 


All  available  recordings  at  the  stations  EKA,  GBA,  WRA  and  YKA  for 
the  events  in  the  NEIS  section  of  table  1  were  phased  to  obtain  the  optimum 
"best  beam".  Because  the  time  code  was  absent  on  WRA  tapes  for  several  months 
from  17  May  1976,  many  of  the  seismograms  could  not  be  found  at  this  station. 
Table  1  shows,  for  each  earthquake,  which  array  seismograms  are  available  and 
whether  or  not  the  P-wave  seismogram  was  detected.  "Detections"  include  any 
cases  where  a  phased  arrival  could  be  seen  on  the  cross-correlogram  of  two  semi¬ 
sums  (filtered  between  1  and  2  Hz)  since  this  has  a  lower  detection  threshold 
than  the  phased  array  sum. 


Examination  of  the  seismograms  revealed  that,  as  expected,  the  three 
events  with  m^  greater  than  6  (events  1,  2  and  32  of  table  1  -  figures  25  and  43) 
show  highly  complex  short  period  seismograms,  and  so  their  mechanisms  are  not 
considered  further  in  this  report.  Examination  of  the  remainder  revealed  that 
most  consist  of  one  or  more  discrete  arrivals  which  could  be  provisionally 
identified  as  P,  pP  and  sP.  (In  most  cases,  the  focal  mechanism  computations  will 
add  strength  to  these  interpretations.)  On  the  basis  of  this,  the  events  were 
classified  into  five  categories  as  follows:- 


Category  1: 


Category  2: 


Category  3: 


Category  4: 


Category  5; 


Earthquakes  with  overloaded  or  extremely  complex  seismo¬ 
grams  (the  three  large  events  already  eliminated). 

Earthquakes  with  clear  discrete  phases  identifiable  on  at  least 
two  (but  generally  three  or  four)  seismograms,  with  consistent 
relative  arrival  times  at  each  station. 

Earthquakes  with  at  least  two  well-recorded  seismograms  but 
with  complex  waveforms  or  inconsistent  phase  arrival  times 
between  stations. 

Earthquakes  with  only  one  seismogram  available  with  good 
signal  to  noise  ratio,  ie,  a  limitation  is  imposed  by  the 
unavailability  of  one  or  more  arrays,  rather  than  by  signal 
amplitude. 

Low  m.  earthquakes  with  low  signal  to  noise  ratio  seismograms, 
or  below  the  detection  threshold. 


Events  are  classified  in  this  way  in  table  1,  and  in  figure  3  which 
shows  the  relationship  between  magnitude,  category  and  time.  This  figure  gives  a 
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preliminary  indication  of  the  magnitude  threshold  above  which  the  relative 
amplitude  method  might  give  useful  results;  this  will  be  discussed  later. 


Our  criteria  for  classifying  these  earthquakes  on  the  basis  of  their 
array  waveform  can  be  seen  from  examination  of  the  seismograms.  Those  for 
Category  2  are  shown  in  figures  5(b)  to  23(b),  and  those  lor  Categories  1,  3,  4  and 
5  are  shown  in  figures  25  to  67,  together  with  correlograms  filtered  between  1 
and  2  Hz.  The  correlograms  indicate  the  arrival  of  phased  energy  which  is  useful 
in  cases  where  the  signal  to  noise  ratio  is  low.  Note  that,  for  events  in  Category 
1,  some  array  seismograms  were  overloaded  and  are  not  shown. 

For  a  given  earthquake  there  are  several  factors  which  might  make  it 
impossible  to  apply  the  relative  amplitude  method,  and  the  purpose  of  the  above 
classification  is  to  establish  the  relative  importance  of  these.  Thus,  the 
Category  1  earthquakes  are  so  large  that  their  short  period  seismograms  are 
highly  complex  and  possibly  overloaded,  with  no  identifiable  phases.  Indeed  we 
expect  the  duration  of  source  radiation  to  be  outside  the  pass  band  of  the 
seismometer  for  m^  greater  than  about  6.0.  Of  the  earthquakes  below  this  limit, 
we  are  prevented  from  applying  the  method  either  because  one  or  more 
seismograms  are  too  complex  (Category  3),  or  because  some  array  stations  were 
not  operating  (Category  4),  or  because  the  earthquakes  were  too  small  (Category 
5).  It  is  important  to  emphasise  that  of  these  categories,  only  Category  3 
represents  a  fundamental  limitation  on  the  applicability  of  the  method,  whereas 
Category  4  requires  only  that  recording  gaps  be  reduced,  and  Category  5  requires 
that  the  signal  to  noise  improvement  be  increased  (ie,  that  the  array 
instrumentation  or  processing  be  improved  to  lower  the  detection  threshold).  It  is 
a  remarkable  feature  of  this  suite  of  data  that  only  four  events  are  in 
Category  3. 

5.  SOURCE  MECHANISM  DETERMINATIONS 


The  relative  amplitude  method  was  applied  to  all  the  earthquakes  in 
Category  2  of  table  1;  the  assumed  phase  identifications  and  amplitude  bounds 
are  included  in  figures  5(b)  to  23(b),  and  figures  80(b)  and  81(b).  We  must  always 
be  aware  of  the  possibility  that  the  direct  P  wave  is  nodal,  in  which  case  care 
must  be  taken  not  to  misidentify  the  second  arrival  (probably  pP)  as  P.  Although 
a  large  apparent  travel  time  residual  of  two  or  three  seconds  would  point  towards 
this  possibility,  examination  of  the  correlogram  is  more  likely  to  reveal  the 
presence  of  a  small  amplitude  earlier  arrival  with  the  correct  azimuth  and  phase 
velocity.  Although  correlograms  and  time  codes  are  not  shown  in  figures  5(b)  to 
23(b),  they  were  in  all  cases  examined  carefully,  and  one  example  of  a  nodal  P 
observation  which  was  detected  unambiguously  from  the  correlogram  is  the  EKA 
observation  of  Event  9  (figure  9(b)),  which  is  discussed  in  reference  (5).  Other 
examples  of  nodal  P  observations  are  the  GBA  record  of  Event  41  (figure  18(b)); 
the  EKA  record  of  Event  8  (figure  8(b)),  and  the  YKA  records  of  Event  8  (figure 
8(b))  and  Event  36  (figure  46).  We  further  note  that  the  misidentification  of 
phases  must  give  rise  to  the  incorrect  specification  of  amplitude  bounds,  with 
consequent  loss  of  compatibility  with  the  double  couple  source  model.  An 
example  of  this  is  the  GBA  record  of  Event  34  (figures  17(b)  and  24(b)  -  see 
later). 


The  width  of  the  bounds  on  each  sP  phase  is  generally  wider  than 
those  of  P  and  pP,  in  order  to  allow  for  several  extra  uncertainties  in  sP 
amplitude  (5).  The  most  important  of  these  is  that  the  surface  reflection 
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coefficient  of  sP  is  highly  sensitive  to  the  near-surface  velocity  structure;  the 
specified  sP  bounds  all  include  extra  width  to  allow  for  this  uncertainty,  as  it  is 
assumed  that  we  have  no  knowledge  of  shallow  velocities. 

The  result  of  applying  the  relative  amplitude  method  to  each  of  these 
earthquakes  is  expressed  graphically  using  a  vectorplot  which,  for  each 
earthquake,  shows  those  source  orientations  which  are  compatible  with  a  given 
set  of  relative  amplitude  observations.  Interpretation  of  a  vectorplot  in  terms  of 
types  of  fault  and  angle  of  slip  is  aided  by  figure  4  (taken  from  reference  (3)). 

Each  vectorplot  presented  here  shows  those  ranges  of  orientations 
which  are  compatible  with  all  the  relative  amplitude  observations  for  one 
earthquake.  An  initial  set  of  vectorplots  was  generated  for  all  the  Category  2 
seismograms  with  the  assumption  that  there  was  no  energy  loss  suffered  by  pP  or 
sP  at  discontinuities  above  the  source,  except  at  the  free  surface.  A  source 
layer  velocity  of  6.1  km/s  was  assumed  for  the  purpose  of  calculating  ray  takeoff 
angles.  These  vectorplots  are  shown  in  figures  5(a)  to  23(a),  and  80(a)  and  81(a), 
and  represent  realistic  results  for  these  shallow  earthquakes  if  nothing  is  known 
about  the  overlying  velocity  structure. 

Several  general  conclusions  can  be  drawn  from  examination  of  these 
vectorplots.  First,  none  of  these  earthquakes  has  failed  to  yield  orientations 
which  are  compatible  with  a  double  couple  source  type.  Second,  the  number  of 
compatible  orientations  generally  reflects  the  level  of  noise  on  the  observed 
seismograms  -  the  better  quality  recordings,  with  larger  m^,  have  orientations 
which  are  more  closely  constrained.  Third,  the  fraction  of  orientations 
compatible  with  each  set  of  seismograms  is  "small"  in  all  cases  except  two, 
which  are  Events  27  and  49.  These  events  give  only  a  poor  constraint  in 
orientation,  and  are  shown  in  figures  80  and  81.  Table  3  shows  the  fraction  of 
orientations  compatible  with  each  set  of  data  -  expressed  as  the  fraction  of 
incompatible  orientations  in  real  space,  or  "significance"  (4).  The  values  indicate 
that  typically  less  than  1%  of  orientations  are  compatible  with  the  three  or  four 
array  seismograms  used. 

Although  the  above  results  exclude  the  effects  of  velocity  structure, 
theoretical  seismograms  can  be  used  to  place  constraints  on  the  possible  velocity 
structures  above  the  source  and  when  an  appropriate  velocity  structure  is 
allowed  for  in  the  relative  amplitude  calculations,  the  surface  reflections  are 
corrected  for  amplitude  loss  upon  passage  through  shallow  discontinuities.  Table 
2  shows  a  realistic  velocity  structure  inferred  partly  from  examination  and 
modelling  of  the  observed  seismogram  (details  are  given  below),  and  the  results 
obtained  with  the  relative  amplitude  program  using  the  same  measurements  as 
above  are  shown  in  figures  5(c)  to  23(c),  and  80(c)  and  81(c),  which  are  shown  in 
addition  to  figures  5(a)  to  23(a),  etc,  for  comparison.  Again  the  fraction  of 
orientations  compatible  with  each  is  shown  in  table  3. 

For  a  given  set  of  data  we  would  expect  the  significance  to  decrease 
(ie,  the  fraction  of  orientations  compatible  to  increase)  as  the  representation  of 
structural  parameters  becomes  nearer  the  true  values.  Pearce  (5)  found  this  to  be 
true  for  an  earthquake  in  the  Gulf  of  Suez.  In  the  present  sequence,  12  out  of  22 
earthquakes  show  a  decrease  in  compatibility  (table  3)  which  may  indicate  that 
the  structure  used  is  inaccurate.  However,  it  is  important  to  emphasise  that  the 
solutions  are  extremely  stable  under  change  of  velocity  structure;  the  fact  that 
similar  solutions  are  obtained  in  both  cases,  and  that  all  the  sets  of  data  show 
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some  compatible  orientations,  indicates  that  this  matter  is  not  crucial  to  the 
identification  of  the  radiation  patterns  as  double  couples. 


All  the  vectorplots  presented  allow  a  range  of  thrust  mechanisms 
with  a  greater  or  lesser  component  of  strike  slip  motion;  in  six  of  the  less  well 
constrained  solutions  a  small  population  of  very  different  orientations  -  either 
near-vertical  strike  slip  or  vertical  dip  slip  -  are  also  allowed.  It  is  important  to 
understand  the  significance  of  these  different  patterns  of  allowable  orientations 
shown  in  each  vectorplot,  and  to  explain  why  these  patterns  occur.  It  is  also 
important  to  interpret  the  angular  bounds  placed  on  the  source  orientations  and 
to  see  in  what  sense  the  orientations  are  well  constrained  and  in  what  sense  they 
are  poorly  constrained.  In  order  to  assist  in  the  qualitative  interpretation  of  the 
vectorplots  it  would  be  helpful  if  each  of  them  was  reproduced  by  means  of  a 
conventional  lower  hemisphere  equal  area  plot,  on  which  the  same  set  of 
acceptable  orientations  could  be  superposed.  Such  a  composite  lower  hemisphere 
equal  area  plot  has  been  produced  which  is  equivalent  to  each  of  the  vectorplots 
in  figures  5(c)  to  23(c)  (ie,  those  with  allowance  for  velocity  structure).  These  are 
shown  in  figures  5(d)  to  23(d).  Because  individual  fault  plane  solutions  cannot 
easily  be  distinguished  when  many  are  superposed,  these  composite  plots  have 
been  produced  at  a  10°  rather  than  a  5a  angular  increment;  this  reduces  the 
number  of  plotted  mechanisms  by  a  factor  of  eight  (but  introduces  the  possibility 
that  small  acceptable  ranges  of  orientations  may  be  missed  altogether).  In  order 
to  enable  each  fault  and  auxiliary  plane  to  be  related,  each  pole  is  indicated  by  a 
small  square. 

Considering  first  the  allowable  thrust-type  orientations,  it  is  seen 
that  some  of  the  earthquakes  (eg,  Event  23  -  figure  12)  give  a  perfect  thrust  with 
no  significant  strike  slip  component  and  a  well-constrained  strike  running  north¬ 
west  to  south-east.  Some  other  earthquakes  (eg,  Event  48  -  figure  19)  have  a 
similar  strike  less  well  constrained,  where  all  orientations  have  a  similar  and 
significant  component  of  strike  slip.  (As  usual,  the  two  distinct  acceptable 
regions  in  figure  19(c)  denote  the  interchange  of  fault  and  auxiliary  planes.)  Still 
other  earthquakes  (eg,  Event  9  -  figure  9)  allow  orientations  with  almost  any 
strike  for  certain  well-defined  slip  angles,  even  though  the  fraction  of 
orientations  which  is  compatible  is  no  larger.  Here  it  is  important  to  realize  that 
it  is  the  "shape"  of  the  acceptable  regions  in  orientation  space  which  has 
changed  -  the  separate  bounds  on  the  slip  angle,  dip  and  strike  are  not  themselves 
indicators  of  the  degree  of  constraint  on  the  orientation,  since  the  acceptable 
range  may  be  of  any  shape.  It  is  seen  in  figure  9(d)  that  the  range  of  acceptable 
orientations  represents  a  fault  plane  which  is  closely  defined  in  dip  and  strike  if 
it  is  assumed  to  be  the  north-west  to  south-east  striking  nodal  plane. 

All  the  above  source  orientations  are  fundamentally  similar  if  we 
assume  the  north-west  to  south-east  striking  nodal  plane  is  the  fault  plane,  but 
there  are  six  earthquakes  which  give  near-vertical  strike  slip  orientations  as  an 
entirely  separate  region  of  acceptable  solutions  in  addition  to  the  thrust  type  - 
on  either  or  both  of  the  vectqrplots.  This  happens  for  Events  3,  8,  19,  22,  26  and 
28  (figures  5,  8,  10,  11,  13  and  14  respectively).  This  is  of  importance  since  these 
two  types  of  orientation  are  likely  to  have  very  different  geophysical 
implications.  A  closer  examination  of  the  vectorplots  shows  that  this  is  just  an 
extreme  case  of  other  earthquakes  where  there  is  an  "annulus"  of  acceptable 
orientations  around  the  central  thrust  orientations,  and  in  some  cases  where  this 
second  annulus  is  widely  outside  the  central  thrust  area,  it  effectively  cuts 
through  the  vertical  strike  slip  or  vertical  dip  slip  parts  of  the  vectorplot. 
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In  order  to  explain  this  we  first  note  that  it  is  the  thrust  type  faults 
which  persist  throughout  all  the  sequence  of  earthquakes,  and  that,  for  most  of 
the  earthquakes,  these  other  fault  types  are  specifically  excluded.  A  limitation 
of  the  present  data  is  the  azimuthal  distribution  of  the  four  arrays  (figure  1(b)) 
which  sample  two  small  ranges  directly  opposite  each  other.  Hence,  if  most  or  all 
four  arrays  show  positive  first  motion  for  P  and  a  negative  for  pP,  the  first 
motion  data  alone  do  not  allow  only  thrust  solutions,  but  vertical  strike  slip 
solutions  can  also  be  oriented  so  that  all  four  arrays  lie  in  positive  quadrants,  and 
three  of  them  can  be  arranged  in  positive  quadrants  for  near-vertical  dip  slip  if 
one  or  two  pP  polarities  are  allowed  to  be  positive.  For  the  larger  earthquakes, 
bounds  on  the  relative  amplitudes  are  easily  able  to  eliminate  these  alternatives, 
but  for  some  of  the  smaller  events,  where  the  relative  amplitude  ranges  are 
wide,  or  where  not  all  first  motions  can  be  read,  then  some  of  these  alternatives 
may  remain. 

The  conclusion  from  this  is  that,  if  only  four  good  stations  are  to  be 
used,  then  it  is  vital  to  secure  a  good  azimuthal  coverage  when  looking  near  the 
detection  threshold,  where  there  may  be  little  relative  amplitude  information 
available  from  each  station. 

There  is  one  event  (Event  34  -  figure  17),  where  the  GBA  seismogram 
cannot  be  interpreted  unambiguously,  and  the  correlogra.n  shows  a  possible  nodal 
P  pulse  before  the  high  amplitude  pulse  which  is  interpreted  as  P  in  figure  17(b). 
This  ambiguity  is  not  resolved  by  examination  of  the  ISC  predicted  arrival  time, 
which  lies  between  the  two  possible  P  waves,  so  a  second  application  of  the 
relative  amplitude  method  was  attempted,  with  the  assumption  that  the  GBA  P 
wave  is  nodal.  This  second  interpretation  also  yielded  compatible  orientations,  as 
shown  in  figure  24,  where  one  orientation  (figure  79)  has  been  chosen  to  generate 
theoretical  seismograms. 

Sykes  and  Burdick  (6)  computed  focal  mechanisms  for  two  of  the 
large  events  (Events  1  and  32  of  table  1)  from  first  motions  and  modelling  of  long 
period  P  waves.  It  is  interesting  to  note  that  their  independently  determined 
solutions  for  these  two  large  earthquakes,  shown  in  figure  87,  give  the  same  type 
of  thrust  mechanism  as  is  obtained  here  for  the  smaller  events  using  the  relative 
amplitude  method. 

6.  THEORETICAL  SEISMOGRAMS 

It  is  noted  that,  within  the  acceptable  range  of  orientations,  all 
orientations  are  equally  probable,  as  all  are  compatible  with  the  data  to  within 
the  limits  imposed  by  noise.  In  view  of  this,  no  attempt  is  made  to  select  one 
orientation  in  preference  to  others.  However,  for  the  purpose  of  computing 
theoretical  seismograms  one  compatible  orientation  is  chosen  from  figures  5(c) 
to  23(c)  for  each  earthquake,  and  the  chosen  orientation  is  in  each  case  arrowed 
on  the  vectorplot.  Theoretical  seismograms  are  generated  using  the  method  of 
Hudson  (7,8)  and  Douglas,  Hudson  and  Blarney  (9),  with  the  velocity  structure  of 
table  2,  and  a  focal  depth  which  agrees  with  the  relative  times  of  P,  pP  and  sP  on 
the  observed  seismograms.  These  theoretical  seismograms  are  shown  in  figures 
5(b)  to  23(b).  All  the  earthquakes  are  assumed  to  be  "small",  in  the  sense  that 
their  pulse  duration  is  less  than  1  s,  and  a  source  radius  of  1.0  km  is  used  for  each 
(table  2.) 


The  value  of  these  theoretical  seismograms  is  that  they  offer  a  means 
of  investigating  how  well  other  features  of  the  observed  seismograms  can  be 
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reproduced,  in  addition  to  the  relative  amplitudes  of  P,  pP  and  sP  which  are 
implicit  in  the  source  orientation.  For  each  station  we  choose  an  anelastic 
attenuation  factor  of  either  t*  =  0.2  s  or  0.6  s  -  whichever  yields  a  dominant 
period  closer  to  that  observed.  (A  more  accurate  matching  of  t*  values  would  not 
be  relevant  to  mechanism  determinations.)  Values  of  t*  are  shown  in  table  3. 
Generally  speaking,  there  is  little  difference  in  the  dominant  periods  either 
between  the  events  or  between  the  stations.  Exceptions  to  this  are  Event  52 
(figure  20(b))  which  has  a  significantly  longer  period  at  YKA  compared  with  EKA, 
and  Event  31  (figure  16(b))  where  an  uncharacteristically  long  dominant  period  is 
observed  at  all  four  stations,  together  with  a  rather  different  pulse  shape.  We 
would  expect  these  effects  to  be  associated  with  larger  sources,  and  indeed  these 
two  events  are  larger  than  most,  with  of  5.3  and  5.1  respectively  (table  2). 

We  note  that,  especially  for  the  seismograms  with  high  signal  to  noise 
ratios,  we  have  succeeded  in  reproducing  the  pulse  shapes  of  P  and  pP.  In 
particular,  the  first  half  cycle  of  pP  is  typically  suppressed;  this  is  most  apparent 
in  the  YKA  seismograms  and  is  reproduced  in  the  theoretical  equivalent. 
Reflection  of  P  energy  from  the  base  of  shallow  sedimentary  layers  gives  rise  to 
small  amplitude  pP  type  arrivals  which  precede  true  pP  by  less  than  a  second, 
thereby  serving  to  increase  the  rise  time  of  the  pP  pulse,  with  the  observed 
effect  on  the  narrow  band  seismogram.  Furthermore,  an  arrival  between  P  and 
pP  is  also  modelled  as  a  reflection  from  a  sub-surface  layer  (see  for  example 
Event  9  (figure  9(b)). 

Another  general  observation  to  be  made  is  that,  although  seismo¬ 
grams  observed  at  the  same  station  for  different  earthquakes  in  this  sequence  in 
many  cases  have  very  different  relative  amplitudes  between  their  phases,  the 
source  orientations  are  all  of  the  same  general  type  -  with  predominantly  dip  slip 
thrust  motion.  Thus,  several  EKA  observations  have  a  nodal  direct  P,  while 
others  show  a  high  amplitude  direct  P  wave.  We  conclude  that  apparently  very 
different  sets  of  seismograms  can  arise  from  very  similar  source  orientations. 
The  reason  for  this  is  that  the  amplitude  of  a  phase  may  become  large  after  only 
a  few  degrees  of  rotation  from  a  nodal  plane. 

Figure  78  shows  a  composite  lower  hemisphere  plot  of  the  source 
orientations  used  in  the  models.  Since  each  was  chosen  independently,  it  gives  a 
general  indication  of  the  strong  tendency  towards  north-west  to  south-east 
striking  thrust  mechanisms,  although  this  tendency  as  shown  in  figure  78  is  only 
qualitative  in  view  of  the  fact  that,  for  each  event,  one  orientation  was  chosen 
from  all  those  which  were  compatible. 

In  the  cases  where  seismograms  with  low  signal  to  noise  ratio  are 
used,  the  confidence  limits  on  the  source  orientation  must  depend  mainly  upon 
the  uncertainty  imposed  upon  the  signal  amplitudes  by  the  noise.  This  uncertainty 
can  be  estimated  by  the  addition  of  noise  with  the  observed  frequency  spectrum 
to  a  theoretical  seismogram,  in  the  way  described  by  Pearce  and  Barley  (10),  who 
generated  whole  series  of  such  noisy  theoretical  seismograms  in  order  to 
establish  the  range  of  effects  that  a  given  character  and  amplitude  of  noise  could 
have  on  the  observed  signal. 

It  is  instructive  to  apply  this  method  in  the  present  situation,  in  order 
to  check  that  the  range  of  effects  obtained  is  compatible  with  the  amplitude 
bounds  which  have  been  placed  on  the  phases.  Figures  68  to  77  show  ten  cases; 
each  shows  the  observed  seismogram,  followed  by  the  theoretical  seismogram 
and  six  randomly  chosen  noisy  theoretical  seismograms. 
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MAGNITUDE  AND  THE  DETECTION  THRESHOLD 


Figure  3  showed  the  relation  between  m.  and  "Category",  when  the 
earthquakes  were  classified  as  in  section  4.  The  successful  application  of  both 
the  relative  amplitude  method  and  modelling  to  the  Category  2  earthquakes  has 
now  confirmed  the  classification  of  these  earthquakes  in  all  but  two  cases 
(Events  27  and  49)  which  did  not  yield  "well  constrained"  source  orientations. 
Neglecting  the  our  Category  3  earthquakes,  we  conclude  that  16  out  of  30 
events  with  m,  between  4.5  and  6.0  have  been  successfully  modelled  as 
earthquakes,  ana  a  further  one  (Event  26)  was  unable  to  be  processed  only  on 
account  of  data  unavailability.  Of  two  events  with  ISC  of  3.6,  both  were 
detected  at  YKA,  and  one  was  detected  at  all  four  arrays  (see  table  1.) 

The  four  Category  3  earthquakes  (ie,  with  complex  seismograms) 
deserve  special  mention.  They  ail  occurred  within  two  days  (see  table  1  and 
figure  3)  which  suggests  that  they  are  in  some  way  related.  Furthermore,  the 
GBA  and  EKA  seismograms  for  these  earthquakes  are  complex,  while  the  YKA 
records  show  the  characteristic  two-pulse  waveform  observed  for  all  the  other 
earthquakes,  albeit  with  some  superposed  complexity.  The  four  GBA  records  have 
similarities  in  their  complexity,  as  do  the  EKA  records.  These  four  earthquakes 
therefore  constitute  a  highly  anomalous  sub-set  of  this  sequence,  but  their 
behaviour  is  unexplained. 

8.  EVENT  LOCATIONS 


In  this  section  we  ask  whether  there  is  any  significant  effect  upon  the 
distribution  or  scatter  of  the  earthquake  sequence  when  the  whole  suite  is 
relocated  using  joint  epicentre  determination  (JED)  (11).  For  this  purpose  the 
depths  of  those  earthquakes  which  are  modelled  are  restrained  to  the  values  used 
in  the  model,  which  are  effectively  determined  from  the  P  to  pP  time.  It  must  be 
remembered,  however,  that  the  identifications  of  pP  have  been  supported  by  the 
success  achieved  in  obtaining  compatible  orientations  when  the  relative 
amplitude  method  was  applied.  In  view  of  this,  the  same  interpretations  were 
placed  upon  those  seismograms  for  earthquakes  which  were  not  modelled,  so  that 
many  more  events  were  given  a  restrained  depth.  In  addition  an  origin  time  and 
hypocentre  were  restrained. 

Figure  82  shows  the  result  including  all  the  ISC  bulletin  data  (with 
Jeffreys  weighting),  with  the  starting  approximation  of  the  NE1S.  Figure  83 
shows  the  same  data  but  with  the  starting  approximation  of  the  ISC.  As  we 
expect,  the  relocated  epicentres  are  the  same  on  each  figure,  to  within  the 
confidence  ellipse  of  the  location;  figure  83  has  some  extra  events  to  the  south 
that  were  not  reported  in  the  NEIS  bulletin. 

It  is  of  interest  to  learn  how  well  these  events  can  be  located  using 
only  the  four  arrays,  which  offer  relatively  accurate  arrival  time  measurements 
from  phased  sums.  Lilwall  (personal  communication)  has  found  that,  for  presumed 
nuclear  explosions  in  East  Kazakhstan,  the  four  arrays  give  as  good  a  location  as 
is  obtainable  using  the  worldwide  data  contained  in  the  ISC  bulletin.  Since  no 
convergence  can  be  obtained  if  the  data  are  underdetermined,  we  can  only  use 
those  earthquakes  for  which  arrival  times  can  be  measured  at  three  or  four 
arrays  -  three  arrival  time  observations  and  a  restrained  depth  provide  four 
measurements  for  the  four  degrees  of  freedom,  namely  latitude,  longitude,  depth 
and  origin  time. 
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Figure  84  shows  these  earthquakes  relocated  using  the  four  arrays 
with  the  NEIS  epicentres  as  starting  approximation;  the  same  epicentre  is 
restrained  as  before.  These  relocated  epicentres  are  compared  in  figure  85  with 
those  obtained  using  the  ISC  bulletin  data  for  the  same  earthquakes  (figure  82). 
In  most  cases  the  array-relocated  epicentres  are  within  the  confidence  limits  of 
those  determined  using  the  ISC  data  (see  table  4),  and  a  major  part  of  the 
difference  is  made  up  of  a  systematic  shift  due  to  the  restrained  epicentre 
(figure  85). 


It  is  interesting  to  note  that  the  P  to  pP  times  give  well  determined 
depths  which  vary  between  5  and  30  km  for  different  events.  This  contrasts 
markedly  with  the  ISC  depths,  many  of  which  have  gross  errors,  especially  when 
the  Quetta  determinations  have  been  used  (see  figure  2).  We  now  ask  whether  the 
relocated  epicentres  exhibit  any  three-dimensional  spatial  trend  when  these 
depths  are  used.  Figure  86  shows  north -south  and  east-west  sections  for  those 
epicentres  computed  in  figure  83  for  which  pP-P  depths  were  used.  No  significant 
spatial  trend  is  apparent. 

9.  CONCLUSIONS 

The  observed  seismograms  presented  in  this  report  are  remarkable  in 
that  the  P  waveforms  comprise  one  or  two  discrete  pulses,  which  have  been 
identified  as  P  and  its  surface  reflections.  Furthermore,  this  behaviour  is 
maintained  for  upwards  of  fifty  earthquakes. 

Application  of  the  relative  amplitude  method  has  confirmed  that  the 
phase  identifications  are  compatible  with  the  double  couple  radiation  pattern  for 
over  half  of  these  earthquakes  between  m.  4.5  and  6.0  (details  in  figure  3),  and 
for  some  of  the  earthquakes  down  to  mb  <P.l.  For  19  of  these  earthquakes  well- 
constrained  focal  mechanism  solutions  have  been  obtained,  indicating  thrust 
faults  with  a  predominance  of  south-west  to  south-east  trending  strikes.  In  a  few 
cases  a  small  number  of  near-vertical  strike  slips  or  vertical  dip  slip  faults  were 
also  possible.  The  thrust  type  orientations  agree  with  the  mechanisms  determined 
for  two  of  the  large  shocks,  which  were  computed  independently  by  a  different 
method  (6). 

The  relative  locations  of  these  earthquakes  have  been  determined 
almost  as  accurately  using  arrival  times  measured  from  phased  sums  at  the  four 
arrays,  as  is  possible  using  ISC  bulletin  data  at  all  stations. 

The  seismograms  presented  here  constitute  a  remarkable  sequence  of 
simple  P-wave  signals  which  contain  much  valuable  source  mechanism  informa¬ 
tion.  Considerable  success  has  been  achieved  in  modelling  these  events  using  an 
earthquake  source.  Further  studies  will  establish  whether  these  array  seismo¬ 
grams  allow  any  of  these  events  to  be  modelled  using  an  explosion  source. 
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TABLE  1 


Uzbekistan  Earthquakes  in  1976  (NEtS  Parameters) 
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Additional 

Events  Reported  to  the 

ISC  ( ISC  Parameters) 

1976 

23  April 
17  May 

17  May 

17  May 

13  June 

19- 55-30.0 
03-21-54.0 
03-57-42.0 
05-53-14.3 

20- 49-34.5 

4.4 

*N  (normal)  indicates  that  depth  was  constrained  to  11.0  km  in  the  NETS  computation. 


**Using  correlogram  if  necessary:- 

/  *  p  wave  detected;  ^  *  doubtful  detection; 

x  •  not  detected;  -  *  station  not  recording  (later  WRA  seismograms  could  not  be  found  because 

of  absence  of  time  code) 


TABLE  2 


Modelling  Parameters 

(a)  Source  region  atructure 

The  S-wave  velocity  in  each  layer  is  assumed  to  be  1//?  of  the 
P-wave  velocity. 


Layer 

V 

p 

Density, 

g/cm* 

Thickness, 

km 

1 

3.0 

2.7 

0.5 

2 

4.6 

2.7 

3.0 

3 

6.1 

2.8 

29.0 

4 

8.2 

3.3 

- 

(b)  Source  mechanism  (Savage  (13)  model  with  rupture  amplitude 
falling  oil  to  edge  of  fault) 

Fault  radius:  1  km  (circular) 


Stress  drop: 
Rupture  velocity: 


100  bars 

0.6  x  S-wave  velocity  in  source  layer 


Focal  Mechanisms  and  Anelaatic  Attenuation  Factors 
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*Mot  modelled 


FIGURE  la  Location  Map  of  Earthquake  Region 


41°N 


63*E 


40°N 


64°E 


YKA 


EKA  €*  312-6* 
V  A=  «-5* 

T\ 

N 

S 

S 


£  =  359  2° 
A  *  77  5° 


V - . - 


\  \  - - 

\  v 
•  GAZU  \ 

\ 

\ 

\ 

\ 

^BA  i  «  151  5° 
A  E  29  3° 

0  10  KILOMETRES 


l  \ 

o-\.. 


A  *  8$  8° 


•  BUKHARA 


KEY 

£  AZIMUTH  OF  ARRAY 

A  EPICENTRAL  DISTANCE  TO  ARRAY 

/  DIRECTION  OF  ARRAY  STATION 
/  FROM  REGION  OF  EARTHQUAKE 


INFORMATION  FROM  REFERENCE  12 

_ , _ „ _  FAULT  PLANE  OF 

APRIL  8  EARTHQUAKE 

_ o _ o _  FAULT  PLANE  OF 

MAY  17  EARTHQUAKE 

.  AFTERSHOCK  AREA 

- INTENSITY  IX  ISOSEISMAL 


FIGURE  1b  Map  showing  Azimuths  and  Epicentral  Distances 
to  the  Four  Array  Stations. Supposed  Fault  Planes 
and  Aftershock  Areas  are  from  Reference  12 
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LATITUDE 


SYMBOL  ISC  DETERMINED  DEPTH 

O  0  -  10  km 

□  11 -40  km 

A  41  - 100  km 

A  >  100km 

X  RESTRAINED  TO  33  km 

NOTES: 

THOSE  WITH  DEPTH  >100  km  ARE  CONSTRAINED  ACCORDING  TO  THE  QUETTA  LOCAL  DETERMINATION 

pP-P  TIMES  GIVE  ALL  OEPTHS  AS  INTRACRUSTAL  (SEE  TEXT  I 

OEPTH  RANGES  ARE  INDICATED  BY  DIFFERENT  SYMBOLS  AS  SHOWN 


FIGURE  2  Plot  of  Uzbekistan  Hypocentres  as  Determined  by  the  ISC 


(NEIS) 


X  CATEGORY 


LARGE  MAGNITUDE  EARTHQUAKES  WITH  HIGHLY  COMPLEX 
ARRAY  SEISMOGRAMS 


o 


CATEGORY  2 


AT  LEAST  TWO  WELL-RECORDED  ARRAY  SEISMOGRAMS 
WITH  CLEAR  DISCRETE  PHASES 


O  CATEGORY  3  AT  LEAST  TWO  WELL- RECORDED  ARRAY  SEIMOGRAMS 
WITH  COMPLEX  OR  INCONSISTENT  WAVEFORMS 


A  CATEGORY  4  ONLY  ONE  GOOO  QUALITY  ARRAY  SEISMOGRAM  AVAILABLE 


•  CATEGORY  5  SMALL  MAGNITUDE  EARTHQUAKES  WITH  DATA  OF  POOR 
SIGNAL  TO  NOISE  RATIO 


FIGURE  3  Relation  beween  ,  Earthquake  Category  and  Time 
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Acceptable  orientations  are  plotted  as  vectors  from  the  Cartesian  point 
defining  \jj  and  6,  in  the  direction  of  the  strike  o.  Lower  hemisphere 
stereographic  projections  indicate  the  type  of  fault  plane  orientation 
represented  by  various  combinations  of  iji  and  6,  and  are  shown  oriented  for 
strike  0  =  360°  (northerly).  In  each  case  the  fault  plane  is  shown  by  o 

thick  line  and  the  auxiliary  plane  by  a  thin  line -  .  Shaded 

quadrants  f  ■  are  negative.  Different  parts  of  the  plot  characterize 
various  fault  types,  and  some  of  these  are  shown.  Where  the  interchange 
of  fault  and  auxiliary  planes  yields  a  different  fault  type,  this  is 
shown  in  square  brackets  [  ]. 


FIGURE  4.  METHOD  OF  REPRESENTING  ACCEPTABLE  FAULT  PLANE  ORIENTATIONS  IN 
TERMS  OF  SLIP  DIRECTION  >h,  DIP  5  AND  STRIKE  a,  AS  DEFINED  IN' 
REFERENCE  (~3) 
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8  April  1976 


03-30-49.3 


(yith  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  2530 


FIGURE  6(c) 
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8  April  1976  06-16-50.0 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  674 


FIGURE  7(a) 


EVENT  NO.  7  (CATEGORY  2) 


P:  5  to  10  *ve 
pPi  4  to  7  -ve 
•Pi  2  to  5  *ve/-v« 


8  April  1976 _ 06-16-50.0 
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FIGURE  7(b) 
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(wit h  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ■  625 
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FIGURE  7(d) 

(equivalent  to  figure  7(c)) 
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8  April  1976  09-10-07.7 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  772 


FIGURE  8(a) 


EVEKT  HO.  8  (CATEGORY  2) 
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(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ”  690 


FIGURE  8(c) 
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DEPTH  DEDUCED  BY  MODELLING: 
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8  April  1976  12-03-41.1 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ”  210 


FIGURE  9(a) 


* 
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TOUT  MO.  9  (CATECOBT  2) 


«  April  1976  12-01-41.1 


(a) 
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PICOM  9(b) 
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DIP  Of  FAULT  PLANE 
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8  April  1976  12-03-41.1 

(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  288 


FIGURE  9(c) 
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FIGURE  9(d) 

(equivalent  to  figure  9(c)) 


EVENT  NO.  19: 

DEPTH  DEDUCED  BY  MODELLING: 
MAGNITUDE  (NEIS) 

SOURCE  ORIENTATION  USED  IN  MODEL 
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17  April  1976  12-11-14.2 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  “  296 


FIGURE  10(a) 


DIP  OF  FAULT  PLANE 
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17  April  1976 


12-11-14.2 


(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  337 


FIGURE  10(c) 
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EVENT  NO.  22: 

DEPTH  DEDUCED  BY  MODELLING: 
MAGNITUDE  ^  (NEIS): 

SOURCE  ORIENTATION  USED  IN  MODEL: 


18  April  1976 
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Slip  Angle  -  60° 
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18  April  1976 


22-37-39.7 


(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ”  341 


FIGURE  11(a) 


8VBWT  <0.  22  (CATEGORY  2) 


18  April  1976  22-37-19.7 
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(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  =  4461 


FIGURE  13(c) 


56 


SCHMID’  EC.IJA. 


FIGURF.  n(d) 

(equivalent  to  figure  life)) 


DIP  OF  FAULT  PLANE 


1 3-57-01.0 


EVENT  NO.  28: 

DEPTH  DEDUCED  BY  MODELLING: 
MAGNITUDE  (NEIS) 

SOURCE  ORIENTATION  USED  IN  MODEL: 


24  April  1976 
9.18  km 
4.7 

Dip  -  130° 

Slip  Angle  -  55° 

Strike  -  260° 


S° 

5°  ♦ 


SLIP  ANGLE  IN  FAULT  PLANE  0 

30°  60°  90° 


t20° 


30° 


60° 


J 

90  V 

A 

4 

♦ 


120°  + 

+■ 

+  ■ 
♦ 


150°  + 


-*• 


180° 


♦  + 
4  4 
4  ♦ 

+  4 

4  + 

♦  + 

4-  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  -4 

+  A 

+ 

4  4 

4  4 

4  4 

A  + 

/ 

>  / 
A  A 

4  4 

4  4 


STRIKE 


4  4 


44444  +  44-f 
444444444 
444444444 
4444444*4 
444444444 
444444444 
444444444 
444444444 
444444444 
444444444 
444444444 
444444444 
444444444 


4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 


1  1  i  1  1  1  1  1  1  1 


♦  •.  +  +  +  +  +  111 

+  +  +  +  +  +  +  +  +  + 

+ +vvvv\\\v 


„**#/»*** V V 

>v**vvvv**r 

^  4  4  4  4  4  4  t  t  t 
+  +/  /44tttl 

*  *  *  r  t  r  t  tit 
+  +  +  +  +  +  tt13 
<*4^/+  +  +  +  +  +  +  + 
4444ttttt1 
+  +  +  /  t  t  r  r  t  r 


i  i  i  i  i  i  i  i  i  i  i 

11111114114 

+  +  +  +  +  +  +  +  +  +  ♦ 
N-VVVVVVV  +  +  + 

**********  1 
\  \  \  \  % 

\  \  *  * 

V>>4*44  4  4^4 
///MAAAAf,* 

1  \  \  \ 

1  \  \\ 

3\\\\*^+  +  +  + 

44444444444 

t  i  t 


150° 


4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

1  1  + 

4  4  4 

4  4  4 

4  4  4 

✓  4 
4 

«**  A  4 
4*  4-  4 
4—4  4 

4—4  4 

4  4  4 

4  4  4 

4  4^  4s. 

4s  *s  -K. 
4  K  4^ 
4  4  4 

4  4  4 


180° 


4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4 

4  4— 


4—  4-  4  4 

4s.  4-  4—  4  4 

4s.  4-  4  4  4 

*s  4  4  4  4 

4  4  4  4  4 

4  4  4  4  4 


4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

♦  V  ♦ 

♦  "V  V 

♦  A  V 

4  4  4 

4  V  * 

4  4V 
4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 
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(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  2526 
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24  April  1976  13-57-01.0 

(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  =  4461 


FIGURE  14(c) 
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(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ■  288 
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Number  of  compatible  orientations  *  184 
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9  May  1976  07-51-16.7 

(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ■  1 


FIGURE  16(c) 
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17  May  1976 


04-53-51.7 


(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  737 


FIGURE  17(c) 
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19  May  1976  15-54-45.6 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  304 


FIGURE  18(a) 
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19  May  1976  15-54-45,6 

(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ■  403 


FIGURE  18(c) 


76 


FIGURE  18(d) 

(equivalent  to  figure  18(c)) 
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28  May  1976  14-05-37.3 

(with  no  allowance  for  velocity  structure) 

Number  of  compatible  orientations  *  143 


FIGURE  19(a) 


EVENT  HO.  48  (CATEGORY  2) 


28  Msy  1976  14-05-37.3 
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Observed  st  YKA 


Computed  for  YKA 
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(d)  Computed  for  GBA 
t*  -  0.6 
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FIGURE  19(b) 
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28  May  1976  14-05-37.3 

(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  m  192 


FIGURE  19(c) 


DIP  OF  FAULT  PLANE 
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20  June  1976  23-33-48.8 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  836 


FIGURE  20(a) 


BVEHT  HO.  52  (CATEGORY  2) 


20  June  1976  23-33-48.8 
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FIGURE  20(b) 
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20  June  1976 


23-33-48.8 


(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  809 


FIGURE  20(c) 
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FIGURE  20(d) 

(equivalent  to  figure  20(c)) 


09-49-32.8 
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23  June  1976  09-49-32.8 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  “  218 


FIGURE  21(a) 


86 


EVEHT  HO.  53  (CATEGORY  2) 


23  June  1976  09-49-32.8 


(a)  Observed  at  YKA 


(b)  Computed  for  YKA 
t*  =  0.6 
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pP:  3.0  to  7.0  +ve/-ve 
sP:  0  Co  4.0  +ve/-ve 
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P:  2.0  to  4.0  +ve/-ve 
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sP:  0  to  2.0  +ve/-ve 


FIGURE  21(b) 
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23  June  1976 _ 09-49-32.8 

(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ■  414 


FIGURE  21(c) 


FIGURE  21(d) 
(equivalent  to  figure  21(c)) 
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DIP  OF  FAULT  PLANE 


23-35-38.0 


EVENT  NO.  54: 

DEPTH  DEDUCED  BY  MODELLING: 
MAGNITUDE  ^  (NEIS): 

SOURCE  ORIENTATION  USED  IN  MODEL: 


8  July  1976 
7.96  ka 
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Dip  -  135° 

Slip  Angle  -  100° 
Strike  -  320° 
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8  July  1976 


23-35-38.0 


(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  ■  97 


FIGURE  22(a) 


EVENT  NO.  54  (CATEGORY  2) 


P:  7  to  10  +ve 
pP:  4  to  7  -ve 
aP:  0  to  2  +ve/-ve 
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FIGURE  22(b) 


DIP  OF  FAULT  PLANE 
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(with  allowance  for  velocity  structure) 
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FIGURE  22(c) 
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EVEHT  MO.  58t 

DEPTH  DEDUCED  BY  MODELLINGS 
MAGNITUDE  mfc  (NE1S): 

SOURCE  ORIENTATION  USED  IN  MODEL: 


22  September  1976  21-49-42.6 
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22  September  1976  21-49-42.6 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientations  m  106 


FIGURE  23(a) 


AD-A089  512  ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT  ALDERMASTON  <EN— ETC  F/6  B/ll 

THE  1976  EARTHQUAKE  SEQUENCE  IN  UZBEKISTAN  -  FOCAL  MECHANISMS  D— ETC(U) 
JUL  80  R  S  PEARCE#  H  BAINBRIDGE#  P  F  KEY 

nun  ACCTCtrn  iUOC.A97/fl  n  nn  j  r  _  nna 


gvsirr  HO.  58  (CATEGORY  2) 


22  Scptcabcr  1976 _ 21-49-42.6 
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FIGURE  23(b) 
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22  September  1976 


21-49-42.6 


(with  allowance  for  velocity  structure) 
Number  of  compatible  orientatioas  »  200 


FIGURE  23(c) 
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FOCAL  PLANES  GA2L1 


22  SEPT  76  2149 


SCHMIDT  EQUAL  AREA  PROJECTION 


FIGURE  23(d) 

(equivalent  to  figure  23(d)) 


DiP  OF  FA'jLT  PLANE 
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DCPTB  DEDUCED  BY  MODELLINGS 
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17  May  1976  04-53-51.7 

(with  no  allowance  for  velocity  structure) 
Number  of  compatible  orientation*  *  1612 


FIGURE  24(a) 
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FIGURE  24(b) 
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17  May  1976  04-53-51.7 

(with  allowance  for  velocity  atructure) 
Ruaber  of  compatible  orientationa  *  1323 


FIGURE  24(c) 
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FIGURE  24(d) 

(equivalent  to  figure  24(c)) 
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FIGURE  27 


EVENT  NO.  10  -  CATEGORY  5 


8  April  1976  22-54-17. B 


Observed 


EKA 


Correlogram 


--  A’.v  -  •  ^ 


Cor re  log ram 


105 


EVENT  NO.  11  -  CATEGORY  5 
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FIGURE  29 
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FIGURE  32 
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FIGURE  33 
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(See  alto  Figure  81) 
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MAGNITUDE  (NEIS): 


23  April  1976  20-55-31.7 
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Although  assigned  to  Category  2  the  source  orientation  is  less  well 
constrained  than  others  in  Category  2,  and  the  earthquake  has  not 
been  modelled  (see  text). 
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(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *  5864 


FIGURE  80(c) 
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(with  allowance  for  velocity  structure) 
Number  of  compatible  orientations  *=  4806 


FIGURE  81(c) 
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FIGURE  83  Joint  Epicentre  Determination  Plot 

using  IsC  Data  as  the  Starting  Approximations 
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FIGURE  86  (a)  Depth  plot  for  Gazli  Sequence  (n/s  section) 
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Nodal  planes:- 

(1)  Strike  -  91°,  Dip  =  43 

(2)  Strike  -  281°,  Dip  -  48 


FIGURE  87(a).  SOURCE  ORIENTATION  OF  8  APRIL  1976  02-40-27.0  EARTHQUAKE  (EVENT  1  OF 
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Quantity 

Unit 

Symbol 

Conversion 

Basic  Units 

Length 

metre 

■ 

1  a  -  3.2808  ft 

1  ft  -  0.3048  a 

Mass 

kilogram 

*«* 

1  kg  -  2.204b  lb 

1  lb  -  0.45359237  kg 

1  ton  *  1016.05  kg 

Derived  Units 

Force 

nevton 

N  •  kg  a/a2 

1  N  •  0.2248  lbf 

1  lbf  •  4.44822  N 

Work,  Energy,  Quantity  of  Heat 

Joule 

J  •  N  «R 

1  J  -  0.737562  ft  lbf 

1  J  -  9.47817  »  lO-1*  Btu 

1  J  -  2.3884b  «  10"*  kcal 

1  ft  lbf  -  1.35582  J 

1  Btu  -  1055.0b  J 

1  kcal  -  418b. 8  J 

Power 

watt 

W  -  J/e 

1  W  -  0.23884b  csl/s 

1  cal/s  -  4 . 18b8  W 

Electric  Charge 

coulomb 

C  •  A  a 

- 

Electric  Potential 

volt 

V  -  W/A  -  J/C 

- 

Electrical  Capacitance 

farad 

F  -  A  s/V  *  C/V 

- 

Electric  Resistance 

ohm 

n  -  v/a 
s  -  i  rr1 

* 

Conductance 

slanen 

- 

Magnetic  Flux 

weber 

Wb  »  V  s 

- 

Magnetic  Flux  Density 

teals 

T  -  Wb/*2 

- 

Inductance 

henry 

H  -  V  s/A  •  Wh/A 

Conplex  Derived  Units 

Angular  Velocity 

radian  per 

second 

rad/s 

m/s2 

1  rsd/s  »  0.159155  rev/s 

1  rev/s  ■  b. 28319  rad/s 

Acceleration 

metre  per  square  i 

second 

1  m/s2  -  3.28084  ft/s2 

1  ft/s2  -  0.3048  m/s2 

rad/s2 

Angular  Acceleration 

radian  per 

square 

second 

- 

Pressure 

newton  per 

square 

metre 

N/m2  •  Pa 

1  N/a2  -  145.038  «  I0*6  lbf/ln? 

1  lbf/ln?  -  ...8947b  *  103  N/*2 

bar 

bar  -  10s  N/m2 

- 

1  In.  Hg  -  3 38b. 39  N/a2 

Torque 

newton  metre 

N  a 

1  N  m  -  0.737562  lbf  ft 

1  lbf  ft  -  1.35582  N  m 

Surface  Tension 

newton  per 

metre 

N/m 

1  N/m  •  0.0685  lbf /ft 

1  lbf /ft  -  14.5939  N/m 

Dynamic  Viscosity 

newton  second  per 

square  metre 

N  a/*2 

1  N  s/m2  •  0.0208854  lbf  s/ft1 

1  lbf  s/ft2  -  47.8803  N  s/«‘ 

Kinematic  Viscosity 

square  metre  per 

second 

m2/s 

1  m2/s  -  10.7630  ft2/a 

1  ft2/s  -  0.0929  m2/s 

Thermal  Conductivity 

watt  per  metre  kelvln 

W/m  K 

- 

Odd  Units* 

Radioactivity 

hecquerel 

Bd 

1  Bq  -  2.7027  «  10_n  Cl 

1  Cl  •  3.700  «  10* 0  Be 

Absorbed  Dose 

gray 

Gy 

1  Gy  •  100  rsd 

1  red  •  o.oi  Gy 

Dose  Equivalent 

s lever t 

Sv 

1  Sv  •  100  res. 

1  rest  •  0.01  Sv 

Exposure 

coulomb  per  kllogrem 

C/kg 

1  C/kg  -  3876  R 

1  R  -  2.58  «  lO-o  C/kg 

Rate  of  Leak  (Vacuum  Systems) 

millibar  Utra  per  second 

mb  1/s 

1  mb  •  0.750062  torr 

1  torr  “  1.33322  mb 

•These  terns  are  recognised  terns  within  the  aetrlc  eve tan. 


